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Abstract We review the current state of empirical knowledge of the total budget 
of baryonic matter in the Universe as observed since the epoch of reionization. Our 
summary examines on three milestone redshifts since the reionization of H in the 
IGM, z = 3, 1, and 0, with emphasis on the endpoints. We review the observational 
techniques used to discover and characterize the phases of baryons. In the spirit of 
the meeting, the level is aimed at a diverse and non-expert audience and additional 
attention is given to describe how space missions expected to launch within the next 
decade will impact this scientific field. 



1 Introduction 

Although baryons are believed to be a minor constituent of the mass-energy budget 
of our universe, they have played a dominant role in astronomy because they are 
the only component that interacts directly and frequently with light. Indeed, much 
of modern astrophysics focuses on the production and destruction of heavenly bod- 
ies comprised of baryons. Current observations of baryons extend from our Solar 
System to the very early universe [z > 7), i.e. spanning over 95% of the lifetime of 
our known universe. In terms of cosmology, the principal areas of baryonic research 
include measuring their total mass density, identifying the various elements and 
phases that comprise them, and resolving their distribution throughout the universe. 
In turn, astronomers are increasingly interested in probing the baryonic processes 
that feed galaxy formation and the feedback processes of galaxies that transform 
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and return baryons to the intergalactic medium (IGM). Our mandate from the or- 
ganizers was to review our knowledge of baryons to a broad astronomical audience 
from the epoch of reionization (z w 6 — 10) to the present day. Granted limited time, 
we focus on our empirical knowledge and limit the discussion of theoretical inquiry. 

In the past two decades, advances in telescopes and instrumentation have signif- 
icantly advanced our understanding the distribution of baryonic matter in the cos- 
mos. Modern surveys of the local universe yield an increasingly complete census of 
galaxies and the large-scale structure in which they reside, the distribution of galaxy 
clusters and mass estimates of the hot gas within them, and a view of the diffuse gas 
that lies between galaxies. In this proceeding we review this work, placing particu- 
lar emphasis on techniques related to observing diffuse baryonic phases. Presently 
these techniques are most efficiently pursued at z ^ 3 and z ~ — 0.5, and observa- 
tional constraints on the properties and distribution of baryons are most precise at 
these epochs. It is somewhat embarrassing that there lies a nearly 10 Gyr gap in our 
knowledge of the majority of baryons between these two epochs. Future missions 
need to address this hole in addition to the primary uncertainties remaining for the 
z and z^3 epochs. 

This proceeding is organized as follows: In §|2]we provide a basic introduction 
to the observational techniques used to trace baryons throughout the universe. The 
current best estimates for the total baryonic matter density is reviewed in § |3] We 
review the observational constraints on the main phases of baryons at z = 3 in § H] 
briefly comment on our general absence of knowledge at z = 1 (§|5]l, and then review 
the z ~ missing baryons problem in §|6] We finish with a brief review of planned, 
proposed, and desired mission and telescopes that would significantly impact this 
science. 



2 A Primer on Tracing Baryons 

A complete baryon census poses a diverse set of observational challenges to astron- 
omy, where we most often study baryons by observing the radiation that they emit, 
absorb or reflect. This includes stars, HII regions, cluster gas, planets, debris disks, 
etc. With these observations, one can assess the luminosity, temperature, chemical 
composition, and size of the object under study. In terms of performing a complete 
census and analysis of baryons in the present and past universe, however, this ap- 
proach is currently limited by several factors. First, there is the simple fact that more 
distant objects are fainter. Although recent advances in telescopes, instrumentation, 
and search strategies now provide samples of galaxies at high redshift, these sam- 
ples are sparse and (at z > 1) are comprised of only the brightest objects. There is 
no a priori reason why the budget of baryons would be dominated by the bright 
luminous objects, so these samples do not necessarily account for a large fraction 
of the budget. Second, it is generally difficult to estimate precisely the mass of an 
emitting object. With stars and galaxies, for example, mass estimates generally rely 
on stellar evolution modeling which is subject to many uncertainties (e.g. the initial 
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mass function) that are difficult to resolve even in the local universe. Third, it turns 
out that only a small fraction of the universe's baryons are in the collapsed, lumi- 
nous objects that are most easily detected. This point is especially true in the young 
universe but it even holds in our modern universe. And finally, although diffuse gas 
emits Une and continuum radiation, it is very difficult to detect this radiation even 
from gas in the nearby universe. 

It is the last reason, in particular, that has driven observers to absorption-Une 
techniques for characterizing the mass density, temperature, and distribution of the 
bulk of the baryons in the universe. The absorption-line experiment is standard ob- 
servational astronomy in reverse: one observes a bright, background source (e.g. a 
quasar or gamma-ray burst) with the aim to study the absence of light due to ab- 
sorption by intervening baryons. These basic analytic techniques w ere d eveloped 



primarily to study the interstellar medium (ISM) of our Galaxy (e.g. 1114711 ') via UV 
spectroscopy of bright O and B stars. Baryons with at least one electron will exhibit 
both resonant line absorption (e.g. Lya) and continuum opacity (e.g. the Lyman 
limit feature), primarily at UV and X-ray frequencies. The principal observable of 
an absorption line is its equivalent width Wi, the fraction of light over a spectral 
interval that is absorbed by the gas. The principal physical parameter is the col- 
umn density A^, the number of atoms per unit are^ along the sightline. This is the 
number density equivalence of a surface density. For weak absorption, the column 
density scales linearly with the equivalent width. For stronger (i.e. saturated) lines, 

depends weakly on and may be difficult to estimate. To assess these issues, it 
is highly desirable to obtain data at spectral resolution and to measure directly the 
optical depth profile of an absorption feature. In this manner, one can integrate the 
column density provided accurate knowledge of the atomic data of the transition. 

Table [Tlpresents the principal absorption line features used to study H, He, and 
metals in diverse phases. Because most of these spectral features have rest wave- 
lengths at ultraviolet or higher energies, one must utilize UV and X-ray telescopes to 
probe the past « 7 Gyr of the universe (i.e. z < 1). These satellite missions have lim- 
ited aperture and (often) low-throughput instrumentation which requires very bright 
background sources. Most of the gas, at least at z > 1, is relatively cool (T ^ lO'^K) 
and the observed absorption lines have widths less than « 40 km s^^ To spectrally 
resolve such features, one is driven toward echelle spectrometers and generally the 
largest telescopes on Earth (or in space). In these respects the z^O absorption-line 
experiment is much more challenging than studies of the high z universe, where sev- 
eral of the key features are redshifted into optical pass-bands and accessible to large 
ground-based telescopes. This has led, in a few ways, to a greater understanding of 
baryons in the z ~ 3 universe than our recent past. 

A tremendous advantage of the absorption-line experiment is that it achieves ex- 
tremely sensitive limits for studying diffuse gas. For example, the signal-to-noise 
and resolution easily afforded by current telescopes and instrumentation allows de- 
tection of the HI Lya transition to column densities A^hi ~ lO'^cm^^. This is ten 



' Although we generally view stars and quasars as point sources, they have finite area. One may 
envision the concept of column density like a core sample of the Earth, where spectra resolve the 
location of gas in velocity not depth. 
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Table 1 Key Baryonic Diagnostic Lines and Features 



Line 



Phase 



(K) 



T 




(A) 



(A) 



Lyman- Werner Molecular gas 10-100 

21cm Atomic gas 100-1000 

Lya Atomic+Ionized gas 100-40000 

Ha Ionized gas 10000^0000 

Lyman limit Ionized gas 10000^0000 

Hell Ionized gas 10000^0000 

CrV Ionized Gas 20000^0000 

OVI Warm/Hot Gas 20000-10*' 

OVII,OVIII Hot Gas 10''-10'^ 

NeVIII Hot Gas 10^ 



21.6,18.9 40 8 200 
775 1550 3100 7750 



-1000 2000 4000 1 
21cm 0.7 Ghz 0.4 Ghz 140 MHz 



1216 2400 4800 1.2 

6560 13000 26000 65000 

912 1800 3600 0.9 

304 450 912 0.2 

1550 3000 6000 1.5 

1030 2000 4000 1 



orders of magnitude lower than the surface density of a typical molecular cloud 
in a star-forming galaxy. With coverage of the full Lyman series, one can measure 
the surface density of atomic hydrogen across these ten orders of magnitude and 
probe the densest material to very diffuse environments. By the same token, one 
can study cold H2 {T < lOOK) molecules through Lyman/Werner absorption bands 
while also probing diffuse, hot {T > lO^K) gas via metal-line transitions of O, N, 
and Ne. One thereby traces structures in the universe with a wide dynamic range of 
sizes: molecular clouds with diameters less than 100 pc in the centers of galaxies to 
'voids' spanning tens of Mpc. 

There are several disadvantages, however, to absorption-line studies: (1) the ma- 
jority of key diagnostics have rest-frame UV or X-ray frequencies. At low z, one 
requires a space mission with large collecting area and efficient spectrometers to 
perform the experiment. This is expensive and often technologically infeasible; (2) 
One requires a bright, background UV/X-ray sourc^ subjecting the final analysis to 
biases in the discovery and selection of these sources; (3) Bright UV/X-ray sources 
are rare and sparsely distributed across the sky. Therefore, multiplexing has limited 
value and observational surveys are relatively expensive. Furthermore, with only a 
limited number of sightlines to probe the universe, the densest (i.e. smallest) struc- 
ture are at best sparsely sampled; (4) The technique is limited to the most distant 
luminous objects detected, currently z w 6; (5) dust intrinsic to the gas under study 
may extinguish the background source and preclude the analysis altogether 

Several key aspects of absorption-line research are illustrated in Figure [T] which 
presents the echelle spectrum of quasar J133941.95H-054822.1 acquired with the 
Magellan/MIKE spectrometer |60. This optical spectrum shows the rest-frame UV 
emission of the zqso = 2.980 quasar. One notes the broad Lya-nNV emission line 
superimposed on a underlying, approximately power-law spectrum. All of the ab- 
sorption features, meanwhile, are associated with gas foreground to the quasar The 
most prominent is the continuous opacity at A < 3600A which marks the Lyman 



^ Traditionally, researchers have used quasars although increasingly transient gamma-ray burst 
afterglows are analyzed (e.g. (2Tlll67ll ). 
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limit of a gas 'cloud' at zlls ~ 2.91 . This Lyman limit system (LLS) must have 
an HI column density A^hi > lO'^cm^^ to exhibit a Lyman limit opacity Xll > 1- 
The thicket of narrow absorption lines at A = 3600 to 4850A is dominated by Lya 
absorption from gas with z < Zqso- These lines are termed the Lya forest and they 
trace the intergalactic medium (IGM) of the high z universe. Their HI column den- 
sities are generally less than lO'^ atoms per cm^. Within the Lya forest one also 
notes a very broad absorption line at A « 4360A whose observed equivalent width 
is ~ 50A. This 'redwood' of the Lya forest is a damped Lya system (DLA); its 
very large HI column density, A^hi > 2 x 10^" cm^^, allow the Lorentzian damping 
wings to be resolved in the Lya transition. This feature marks the presence of a 
foreground galaxy and has a larger HI column density than total of all other lines 
in the Lya forest. Finally, one observes a more sparse set of absorption features at 
A > 4900A which are metal-line transitions from ions and atoms of O", Si^, C^, 
C+^, Fe+, etc. These lines enable studies of the metal abundance of our universe 
1 1 13i 135 1 and also probe gas phases that are difficult to examine with Lya alone. 




Fig. 1 Magellan/MIBCE echelle spectrum of the zqso = 2.980 quasar J133941.95+054822.1. The 
data reveal the broad Lya+NV emission line of the quasar imposed on its otherwise nearly power- 
law spectrum. Imprinted on this spectrum is a plethora of absorption features due to gas foreground 
to the quasar. In particular, one identifies continuum opacity at A < 3600A due to a Lyman limit 
system at zlls = 2.97 and a thick of absorption features related to the Lya forest. We also mark 
the damped Lya profile at A r; 4360A of a foreground galaxy {zdla = 2.59). 
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3 The Baryonic Mass Density 



Our first problem is to assess the total budget of baryonic matter in the Universe, 
preferably by independent measurement. The baryonic mass density pi,, generally 
defined relative to the critical density = 9.20 x lO^^^/iyQ gcm^^, Qi, = Pb/Pc 
and /jvo, the Hubble constant in units of 70 km s^' Mpc^\ are fundamental cos- 
mological parameters. The first successful method for estimating 12^ is to compare 
measurements of isotopes for the light elements (H, D, ^He, "^He, Li) against pre- 
dictions from Big Bang Nucleosynthesis (BBN) theory. The relative abundances of 
these isotopes are sensitive to the entropy of the universe during BBN, i.e. the ra- 
tio of photons to baryons. At first, observers focused on He even though Li and D 
are much more sensitive 'baryometers'. This was because He/H was inferred from 
nearby emission-line galaxies where one could obtain exquisite S/N observations. 
It is now realized, however, that the precision of the inferred He/H ratio is limited 
by systemati c un certainties in the modeling of HII regions and the precision of the 
atomic data llOl | . 

The tightest constraints on Qt, using BBN theory come from direct measurements 
of the D/H ratio. Ideally this ratio is measured in a metal-poor gas to minimize the 
likelihood that D has been astrated by stars. The first detectio ns of D , however, were 
from absorption-line studies in the metal-rich Galactic ISM il26ll . These observa- 
tions establish a lower limit to the primordial D/H value of w 1 x 10^^ (e.g. ll89ll ). 
In the past decade, echelle spectrometers on lOm-class telescopes have extended 
the experiment to high redshift lfl6ll . Researchers have searched for the very rare 
'clouds' which are metal-poor and therefore minimally astrated, have large HI col- 
umn density, and are sufficiently quiescent that the higher order lines DI and HI 
Lyman series are resolved (the velocity separation is 82km s^'). Figure |2] shows a 
recent example 19711 . The first few transitions of the HI Lyman series are broadened 
by the damping wings of their 'natural' Lorentzian profiles and a fit to the profile 
establishes A^hi = 10^*'-^^^'' ''^cm^^. The DI transitions are lost within the cores of 
these strong absorption lines but higher order lines show resolved and unsaturated 
DI absorption. Analysis of these lines provide a precise estimation of the DI column 
density. 

The D/H ratio is estimated from the A^(DI)/A^(HI) ratio under the reasonable 
assumption that these isotopes have identical ion configuration^ and differential 
depletion. Allowing for these assumptions, the estimation of D/H is direct: abun- 
dances measurements from absorption-line studies are simply a matter of counting 
atoms. Unlike most other methods to measure abundances, the results are essen- 
tially independent of the local physical conditions in the gas. Figure [3] presents the 
'gold-standard' set of D/H measurements from quasar absorption line surveys of 
high z, metal-poor absorbers. The observations are scattered about log D/H — —4.55 
with a boot-strap dispersion of 0.04 dex. Adopting standard BBN theory, this con- 
strains 7] = ny/ny = 5.8 ±0.7 x 10^'" where iij, and iiy are the number densities of 
baryons and photons respectively. Direct measurements of the cosmic microwave 



' The gas is self-shielding and likely has a high neutral fraction. 
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Fig. 2 H 1 and D 1 Lyman series absorption in the z = 2.70262 DLA towards SDSS 1558-0031 
19711 . For the Lya and Ly/3 transitions, the data is unnormalized and the dashed line traces the esti- 
mate for the local continuum level. The remaining Lyman series transitions are shown continuum 
normalized. The solid green line shows the best single-component fit to the D I and H 1 absorption. 
The estimate of the HI column density comes from analysis of the damping wings present in the 
Lyman a-S lines whereas constraints on the Dl column density come from the unsaturated D I 
Lyman- II transition (~918 A rest wavelength). 



background give ;iy, and one derives a baryonic matter density relative to the critical 
density p,., Qb = Pb/Pc = (0.043 ± 0.003)/i7ol It is a triumph of the BBN theory 
and observational effort that this value has been confirmed by observations of the 
CMB power spectrum ll37[|93ll . 

Nevertheless, a few outstanding issues remain. The dispersion in D/H at high 
z exceed^ that predicted from the reported errors on the individual measurements 
Jt^I- Although this most likely reflects over-optimistic error estimates by the ob- 



Interestingly, the same is true of the Galactic ISM measurements | 63!i|. This result is unli kely to 
be due to observational error and the nature of the dispersion remains an open question I83llll6ll . 
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Fig. 3 Values of the D/H ratio vs. Wm lOTh . The horizontal lines represent the weighted mean and 
jackknife errors of the 6 measurements, logD/H = —4.55 ±0.04. The right hand axis shows how 
the values of D/H translate into values for Qi,h^ using BBN. 



servers, it is worth further exploring models of inhomogeneous BBN 116011 and sys- 
tematic error associated with the astrophysics of H and D ll36ll . The other outstand- 
ing issue concerns an inconsistency between the B BN-predicted Li abundance based 
on the rj value from D/H and observation 1^ 146 1. It is now generally expected that 
this observed Li underabundance will be explained by non-primordial stellar astro- 
physics (e.g. ^^), but a comprehensive model of these processes has not yet been 
developed. 



4 Baryons at z ~ 3 

This section reviews constraints on the baryon census during the few Gyr follow- 
ing reionization, starting with the densest phases and proceeding to the least dense. 
Perhaps not coincidentally (for hierarchical cosmology), this organization also pro- 
ceeds roughly from the lowest cosmological mass density to the largest. Along the 
way, we highlight the key uncertainties in these values and comment on the impact 
of future missions. 
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The majority (but not all) of galaxy samples identified at z > 2 are discovered 
through spectral or photometric signatures associated with ongoing star-formation. 
These include the Lyman break galaxies (see Shapley's presentation), the Lya emit- 
ters (e.g. [50]), sub-mm galaxies lfl9ll . and gamma-ray burst host galaxies (e.g. H^l)- 
These populations dominate star formation activity during the young universe. It is 
likely that they also dominate the stellar mass density in this epoch because 'red 
and dead' galaxies are very rare at this time (e.g. 17711 ). It is difficult, however, to 
estimate the stellar mass density of these star-forming galaxies. This is especially 
true when one is limited to optical and near-IR photometry which corresponds to 
rest-frame UV and blue light, i.e. photons from massive star formation that provide 
little, if any, information on the low-mass end of the IMF. 

To infer the stellar mass, one must adopt a mass to light ratio (M^/L). Standard 
practice is to model the stellar population that gives rise to the photometry and spec- 
troscopy of the galaxy. This modeling is sensitive to the presumed star-formation 
history, reddening by dust intrinsic to the galaxy, and uncertainties in the initial mass 
function. Altogether, these lead to large systematic uncertainties in the stellar mass 
estimate. These issues aside, current surveys do not probe the fainter galaxies of 
what appears to be a very steep luminosity function (e.g. Il2l|]). The contribution by 
Dickinson provides a deeper discussion of these issues. Current estimates of the stel- 
lar mass density at z = 3 range fro m = 2 to 4 x 10^^ hjoMa; Mpc^^ with a sys- 
tematic uncertainty of « 50% ll44lll28ll . This corresponds to 0.005 ± 0.002i2f,/!7o. 
The launch of new satellites with greater sensitivity in near and mid-IR pass-bands 
(Herschel, JWST) will reduce the uncertainties by observing the stellar component 
at redder wavelengths and also by extending surveys to much fainter levels. Even 
with the large current uncertainty, however, it is evident that stars make a nearly 
negligible contribution to the total baryon budget at these redshifts. 



4.2 Molecular Gas 

In the local universe, star formation is associated with molecular gas in the form 
of molecular clouds. Although the link between molecular gas and star formation 
has not been extensively established at z ~ 3, it is reasonable to expect that the 
large SFR density observed at this epoch implies a large reservoir of molecular gas 



within high z galaxies (e.g. II57I1 ). Indeed, m olecular gas is detected in CO emission 



from a number of z > 2 galaxies lll5Ull68ll with faciUties such as CARMA, IRAM, 



and CSO. The observations focus on CO and other trace molecules because the H2 
molecule has no dipole moment and therefore very weak emission lines. Unfortu- 
nately, current millimeter-wave observatories do not afford extensive surveys for 
CO gas in high z galaxies. The Large Millimeter Telescope (LMT) and, ultimately, 
the Atacama Large Millimeter Array (ALMA) will drive this field wide open. These 
facilities will determine the CO masses for a large and more representative sample 
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of high z galaxies and (for a subset) map out the spatial distribution and kinematics 
of this gas. Similar to the local universe, to infer the total molecular gas mass from 
these CO measurements one must apply a CO to H2 conversion factor, the so-called 
'X-factor' . The X-factor is rather poorly constrained even in the local universe and 
nearly unconstrained at high z- Therefore, a precise estimate of the molecular com- 
ponent at high z will remain a challenging venture. 

Although observations of H2 in emission are possible at cosmological redshiftjs^ , 
these measurements are not yet common enough to assess the total molecular gas 
budget at z ~ 3. The H2 molecule does exhibit Lyman/Werner absorption bands 
in the rest-frame UV that are strong enough to be detectable in many astrophys- 
ical environments with UV and optical spectrographs. These ro-vibrational tran- 
sitions have large oscillator strengths and permit very sensitive searches for the 
H? molecule, i.e. H2 surface densities of less than lO'^ particles per cm". To date 
the search for H2 has been targeted toward gas with a known, large HI 
column density (A^hi > lO'^cm^^), i.e. the Lyman limit and damped Lya systems. 
This is a sensible starting point; any sightline intersecting molecular gas in a high 
z galaxy would also be likely to intersect the ambient ISM of that galaxy and ex- 
hibit a large HI column density. Although present surveys are small, it is evident 
that the molecular fraction in these HI absorbers is very low. The frequency of pos- 
itive detections is k, 15% l95ll and the majority of these have molecular fractions 
/(H2) < 10^^ 18 ill . Only a handful of these galaxies show molecular fractions ap- 
proaching the typical Galactic values of even diffuse molecular clouds, /(H2) ^0.1, 
along sightlines with comparable HI column density [26]. 

In hindsight the paucity of H2 detections, especially ones with large /(H2), is 
not surprising. The formation of H2 clouds represents a transition from the atomic 
phase, one that leads to gas with much lower T and much higher density. Conse- 
quently, cold molecular clouds phase exhibit much lower cross-section than when 
the gas is atomic. In turn, these clouds will have significantly lower probability of 
being detected by absorption line surveys. Analysis of the distribution of molecular 
gas in the local universe suggests that fewer than 1 in 1000 quasar sightlines should 
pierce a cold, dense molecular cloud lll78ll . In addition to this probalistic argument, 
these clouds generally have very high dust-to-gas ratios and correspondingly large 
extinction. Therefore, even in the rare cases that a quasar lies behind the molecular 
cloud of a high z galaxy, it may be reddened and extinguished out of the optical, 
magnitude-limited surveys that dominate quasar discovery. It is somewhat curious, 
however, that GRB afterglow spectra also show very low molecular fractions in gas 
that is believed to lie very near star-forming regions II6I I. 

In summary, there are very limited current constraints on the mass density and 
distribution of molecules in the young universe. Furthermore, future studies will fo- 
cus primarily on tracers of H2 and will also be challenged to precisely assess this 
very important phase. At present, our best estimate for the molecular mass density 
may be inferred from the atomic and stellar phases. Allowing for an order of mag- 
nitude uncertainty, we can estimate Q^^^ = (0.001 — 0.1)12^,. Therefore, we suggest 
that the molecular phase contributes less than 10% of the mass density budget at 
z~3. 



Baryons: What, When and Where? 

4.3 Neutral Hydrogen Gas 



11 



In the local universe, galaxies have roughly equal measures of molecular and atomic 
hydrogen gas. The latter phase is primarily studied through 21 -cm emission line sur- 
veys using large radio telescopes. Unfortunately, these radio telescopes do not have 
enough collecting area to extend 21 -cm emission-line observations beyond z ~ 0.2 
ifrgll . Although one can trace the 21 -cm line in absorption to much greater redshifts 
(e.g. lfl4ll66ll ). detections are exceedingly rare. This is due in part to the paucity of 
bright background radio sources and in part because galaxies have a relatively low 
cross-section of cold, neutral gas. For these reasons, neutral hydrogen gas at high z is 
principally surveyed via Lya absorption. A total HI column density of w 10^° cm^-^ 
is required for self-shielding to maintain neutral gas with a standard gas- to-dust ra- 
tio against the extragalactic UV background and local UV sources 116511 . At these 
column densities, the damping wings of the Lya transition are resolved even with 
moderate spectral resolution (FWHM w 2A; Figure [U. This spectral feature lends 
the name 'damped Lya system' to the absorbers with A^hi > 2 x 10^'^cm^^. 
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Fig. 4 The HI frequency distribution f{Niii,X) fo r the damped Lya systems identified in the 
SDSS-DR5 quasar database (mean redshift z = 3.1, ll 14ll ). Overplotted on the discrete evaluation 
of f{Nni,X) is the best fitting F-function which is a reasonable (albeit unphysically motivated) 
description of the observations. 



Wolfe and collaborators pioneered damped Lya research be ginning w ith the Kast 
spectrometer on the 3m Shane telescope at Lick Observatory Il73lll74ll . These op- 
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tical observations cover redshifts z > 1.7 which shifts Lya above the atmospheric 
cutoff. The most recent DLA surveys leverage the tremendous pubhc database 
of quasar spectroscopy from the Sloan Digital Sky Survey (SDSS; Hll). Figured 
presents the frequency distribution f{Niii,X) of A^hi values per liA^ni interval per 
absorption pathlength dX for 719 damped Lya systems ide ntified in the SDSS-DR5 
database (see http://www.ucolick.org/'^xavier/SDSSDLA; ll 1411 '). This distribution 
function, akin to the luminosity function of galaxies, is reasonably well described 
by a Gamma function with a 'faint-end slope' a — —2 and a characteristic column 
density A^* w 10^''^ cm^^. The first moment of f{Nni,X) yields the mass density of 

neutral gas Q,„,„{X)dX = ^^^^ / A^hi/(A^hi,^)^^^ where ^ is the mean molec- 

' N,„i„ 

ular mass of the gas (taken to be 1.3), A',,,,,, is the column density which marks the 
transition from neutral to ionized gas (taken here as 10^'^'^ cm^^), Hq is Hubble's 
constant, and pc is the critical mass density. The faint-end slope of f{NHi,X) is 
logarithmically divergent, but the break at A^* gives a finite density. 

Current estimates for Q„eiit from the SDSS-DR5 database give Q„eut = 0.016 ± 
0.002X2/,. There is a modest decline by a factor of « 2 from z = 4 to z = 2 where the 
latter value is (remarkably) in agreement with the neutral gas mass density observed 
today (§ 16.3b . One may speculate, therefore, that the n eutra l gas mass density has 
been nearly invariant for the past « lOGyr (but see 111 1811 ). This mild evolution 
aside, we conclude that the neutral, atomic gas phase at z = 3 comprises only ~ 
1 — 2% percent of the baryon budget at high z. Cosmological simulations suggest 
that this gas comprises the ISM of galaxies in formation, an assertion supported 
by measurements of DLA clustering with br ight star-forming galaxies 111 ll 12411 and 
the detection of heavy elements in all DLAs 111 1311 . We can conclude that the sum of 
galactic baryonic components is a minor fraction (< 10%) of the total baryonic mass 
density at z ~ 3. This is, of course, a natural consequence of hierarchical cosmology 
(e.g. (921). 

With the DLA samples provided by the SDSS survey, the statistical error on £2„eut 
approaches 10% in redshift bins Az ~ 0.5. At this level, systematic errors become 
important, especially biases associated with the selection criteria and magnitude 
limit of optically bright quasars. The bias which has received the most attention is 
dust obscuration 139119911 . i.e. dust within DLAs redden and extinguish background 
quasars removing them from opt i cal survey s. The dust-to-gas ratio in the ISM of 
these high z galaxies is small 191 16611 . however, and the predicted impact on 
^neur is Only of the order of 10%. This expectation has been tested through DLA 
surveys toward radio-selected quasars ||38, 64], which show similar results, albeit 
with poorer statistical power 

Although we would like to extend estimates of Q„eut to z > 5, the experiment 
is challenged by the paucity of bright quasars (although this may be overcome with 
ORB afterglows) and because the collective, blended opacity of gas outside of galax- 
ies (i.e., the IGM) begins to mimic damped Ly a ab sorption. By z = 6, the mean 
opacity of the universe approaches unity at Lya |l7ll l and the experiment is entirely 
compromised. Of course, current expectation is that the universe is predominantly 
neutral not long before z = 6, i.e. Qneut = ^b- These expectations and the evolution 
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of neutral gas at z > 6 will be tested through 2 1 cm studies of the young universe (see 
contributions by Furlanetto and Loeb for discussions of this reionization epoch). 



4.4 Ionized Gas 

As Figure |6] indicates, the mass densities of stars, molecular gas, and neutral gas 
are unlikely to contribute significantly to 12^ at z ~ 3. Unless one invokes an exotic 
form of dense baryonic matter (e.g. compact objects), the remainder of baryons in 
the e arly universe must lie in a diffuse component outside of the ISM of galaxies 



1 10811 . Indeed, the absence of a complete Gunn-Peterson trough in z < 6 quasars 
demonstrates that the majority of baryons are highly ionized at these redshifts[4o|. 
Because it is impossible to directly trace H+ for the vast majority of the mass in the 
diffuse IGM, we must probe this phase through the remaining trace amounts of HI 
gas or, in principle, Hel and Hell Unfortunately, the resonance-line transitions of He 
have wavelengths that preclude easy detection (Tabl e [H t hough the Gunn-Peterson 



effect has been observed in Hell at z > 2.1, see II781I141I1 ). Therefore, the majority 
of research has focused on the frequency distribution of HI column densities for gas 
associated with the IGM (a.k.a. the Lya forest). 

The modern description for the Lya forest is an undulating density fields with mi- 
nor, but important peculiar velocities and temperature variations. In this paradigm, 
it may be more fruitful to assess distribution functions of the smoothly varying op- 
tical depth of the IGM ||3> [tiI [s^l and compare these against theoretical models 
i Toi 881. Nevertheless, the classical approach of discretizing the Lya forest into in- 



dividual absorption lines does yield valuable insight into the distribution and mass 
density of baryons in the intergalactic medium. This analysis is generally summa- 
rized through the distribution function of HI column densities, f{Nwi,X). Figure |5] 
summarizes our current knowledge of the f{N\{i,X) distribution for both the ionized 
and neutral gas. Quasar absorption line analysis to date has provided an assessment 
of /{Nhi^X) from A^hi = lO'^ to 10^^ cm^^ with an important gap spanning from 
A^Hi ~ 10'^^''^cm^^. At column densities larger than 10"^ cm^", the systems are 
too rare to be measured with existing quasar samples. Detections below lO'^cm^^ 
are limited by S/N and by the difficulty of distinguishing absorption from single 
'clouds' from unresolved blends of weaker and stronger lines (e.g. I^sl), while at 
10^' cm^^, the clouds are unphysically large for typical IGM density and ionization 
conditions. 

Let us focus first on the systems exhibiting the largest A^hi values, the super Ly- 
man limit systems (SLLS; also referred to as sub-DLAs). These absorbers are de- 
fined to have lO'^cm^^ <A^Hi < 10^*^^ cm^^. Similar to the DLAs, one can measure 
the HI column densities of SLLS through fits to the dam ping wing s of the Lya pro- 



rngs 

3i. 



file although higher resolution spectroscopy is required 11981110411 . Current surveys 
indicate that /(A^hi,^) flattens below A^hi ~ 10^° cm"^ with J log //(i log A^ w -1.4 
ioill . If the SLLS are predominantly neutral, then they would contribute only a few 
percent to Queut and a negligible fraction of £2/, (Figure|6]). Most of these absorbers. 
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Fig. 5 The f(NHi,X) distribution function of ionized Nm < 10^" cm"- and neutral A^hi > 
10^" cm"^ at z = 3. The data points at A'hi > 2 x 1 0^° c m~^ correspond to the damped Lya sys- 
tems and are drawn from the SDSS-DR5 analysis of |114!]. The data at 19 < log(A'Hi/ cm^^) < 20.3 
correspond to the super Lyman limit systems (SLLS) and were drawn from the analysis of l93l . 
The dashed line traces their best guess at the frequency distribution for the Lyman limit systems 
with 16 < log(A^Hi/ cm"^) < 19. For the Ly « for est {Nui < lO" cm"-) we adopt the results from 
I73I1 normalized to our assumed cosmology Il45ll . Finally, the dotted line is a spline fit to the Lya 
forest data and the functional form for the LLS given by t98ll . Although the turnover in f{Niii,X) 
lO'"''^ cm"^ may be partially due to incompleteness, it is likely a real effect (see also 

TJ). 



however, are partially ionized with mean ionization fraction < x >w 0.9 1 105 . 1 lOll . 
In Figure|6]we display the baryonic mass density of the SLLS inferred from the HI 
frequency distribution and assuming an A^hi -weighted, average ionization fraction 
< X >= 0.9. We infer that the SLLS contribute on the order of 2% of the baryonic 
mass density at z = 3. The majority of the baryons must reside in yet more diffuse 
and highly ionized gas. 

At column densities A^hi = lO'^ to lO'^cm"^ (the Lyman limit systems; LLS), 
the majority of the Lyman series lines are saturated and the damping wings of Lya 
are too weak to measure. For these reasons, it is difficult to estimate A^hi values 
in this interval. In principle, one can precisely measure A^hi values for the partial 
Lyman limit absorbers with A^hi < 10^^ cm"^ whose opacity at the Lyman limit is 
T912 « 1. To date, howe ver, o nly small samples from heterogeneous surveys exist in 
this range of A^hi (e.g. llOST ). Because the Lyman limit feature is easily identified, 
one more easily derives an integral constraint on f{Niii,X) by surveying the inci- 
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Fig. 6 Estimates of the baryonic mass density (relative to Sli,) for various phases of baryons in the 
z = 3 universe. O n the right side of the figure we show the stellar mass density estimate from several 
recent works l44Lll28ll . The uncertainty in this measurement is dominated by uncertainties in stellar 
population modeling and the presumed IMF. The phases corresponding to logA^ni < 22 refer to 
absorption line systems traced by HI gas. The molecular gas content is very crudely estimated to 
lie near the stellar and atomic mass densities. The neutral at omic phase is traced by the DLAs 
where one estimates a mass density on the order of 1 to 2%i2fc ll 14ll . Gas with Whi < 10^" cm^^ is 
predominantly ionized and one must make significant ionization corrections to estimate the mass 
density from the observed HI atoms (9^. The red data points show the mass density from HI atoms 
alone for the LLS whereas the blue points represent estimates of the total mass density (see the text 
for details). The dashed symbols have weak empirical constraints and the dotted symbols have very 
weak empirical constraints. 



dence of LLS in quasar spectra i{X) = J f{Nm,^)dNm (e.g. lfl49ll ). This integral 
constraint includes gas with A^hi > lO'^cm^^, but is dominated by a bsorbers with 
low A'hi values. Current estimates indicate t{X) =7.5 at z = 3 lll06ll . Analysis of 
the SDSS database will greatly improve this estimate forz > 3.2 (O'Meara et al., in 
prep). Furthermore, complete analysis of the Lya series for LLS will yield tighter 
constraints on/(A^Hi,-'^) within the interval A^hi = lO'^ to 10'''cm^^ (Prochteret al., 
in prep.). |l98(l combined the current observational constraints to predict f{Nni,X) 
across the LLS regime. These are presented as the dashed curve in Figure|5]and the 
uncertainty in this estimation is « 50%. 

Not only is f{N}ii,X) very uncertain for the LLS, so too is the ionization state of 
this gas. For a few systems, analysis of metal-Une transitions hav e con firmed the gas 
is predominantly ionized with ionization fractions x > 0.95 (e.g. Ill 1 ill ). The samples 
are too small, however, to reveal the distribution or trend of ionization fractions with 



16 



Jason X. Prochaska and Jason Tumlinson 



Niii. Unfortunately, there is little insight from theory. Previous work with cosmolog- 
ical simulations has severely underestimated the incidence of Lyman limit systems 
(e.g. lH^) and it is likely that a careful treatment of radiative transfer is required 
to assess these systems that are near the transition from neutral to ionized 17511 . We 
can set a very conservative lower limit to their contribution by taking the f{N}ii,X) 
function of O'Meara et al. and assume x = 0, Qlls > 0. 000412/,. Ionization correc- 
tions will undoubtedly increase this limit by one or more orders of magnitude. Even 
in the case of a 100 times increase, it is very unlikely that LLS dominate the bary- 
onic budget at z 3. Figure |6]presents an estimate of the baryonic mass density of 
the LLS where we have adopted the HI frequency distribution presented in Figure|5] 
and estimated the ionization correction from simple Cloudy calculations assuming 
an io nizat ion parameter \ogU = —2.5, inferred from metal-line analysis of several 



LLS lllllll . This implies a mass density estimate of « 0.05X2/,. 

For column densities A^hi ~ 10^^ to 10''*'^cm^^, where Ly/3 and/or Lya are op- 
tically thin, the frequency distribution is reasonably well characterized by a power- 



law f{Nni,X) oc N^i '^ (e.g. l72l |73[1 ). Observers do report, however, a departure 



from this power-law at A^hi < lO'^cm^^ which is likely not simply a consequence 
of incompleteness in the analysis JtsI ). Although the Lya forest absorbers dominate 
the spectrum of any quasar by number, the total number of HI atoms that they con- 
tribute is negligible (« 4 x 10^^i2/,). From cosmological simulations of the high-z 
universe, we expect that the Lya forest arises in undulating density fields of highly 
ionized gas filling large regions («Mpc) with modest overdensity 5p/p ~ 3. These 
same simulations argue that the majority of the baryonic mass density resides in 
the Lya forest, with a contribution possibly exceeding 0.9 fit, (e.g. lfl5[[88i[l"25l ). 
Empirically, the gas is observed to only exhibit high-ion states of elements like C, 
Si, and O which demonstrate a high degree of ionization. 

Because this gas is optically thin, if one can estimate its volume density and the 
intensity of the ionizing radiation field it is relatively straightforward to estimate 
its io nization fraction and thereby the total baryonic mass density of the Lya forest. 



1 13311 has presented a simple but intuitive prescription for modeling the Lya forest as 
gravitationally bound gas clouds whose size is of order the Jeans length. Within this 
prescription, which matches scaling laws derived from cosmological simulations, 
the ionization fraction x of the gas can be related to the HI column density: 

x = 1 - 4.6 X io-« ( ^ — ^ Tr^-'r-'" (1) 

V2 X 10i3cm-2 y 4 12 w 

where T/s, is the gas temperature in units of 10^ K and Fxi is the photoionization rate 
in units of 10^'^ s^'. 

Photoionization should maintain the gas at a temperature near 2 x lO^K, as 
predicted by the cosmological simulations. Within this prescription, therefore, the 
only major uncertainty is the ionization rate of the extragalactic UV background 
(EUVB) radiation field, Feuvb- At z = 3, the EUVB will have contributions from 
both quasars and star-forming galaxies. The contribution of the former is directly 
measured from the luminosity function of z = 3 quasars. The contribution of galax- 
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ies, however, is far more uncertain. Alth ough recent studies offer relatively well con- 
strained UV luminosity functions (e.g. 112111 ). these are measured at A « 1500A in 
the rest frame, and so must be extrapolated to estimate the flux of ionizing photons. 
Furthermore, one must adopt an escape fraction fesc to estimate the flux emitted 
by these star-f orming galaxies. Current empirical constraints are extremely limited 
l 22Lll38lll42ll : the estimates range from fesc = to w 10%. Because the comoving 
number density of star-forming galaxies greatly exceeds that of quasars, it is pos- 
sible that galaxies contribute as much or more of the ionizing flux to the EUVB. 
Thus this ionizing background is uncertain at high redshift, with a corresponding 
uncertainty in the ionization correction for diffuse ionized gas. 

The most conservative approach toward estimating the baryonic mass den- 
sity of the Lya forest is to assume Feuvb = ^qso- 15811 provide an estimate of 
rgso = 4.6 X 10^'^ s^' at z = 3 with an « 30% uncertainty. Convolving Equation[T] 
(assuming F12 = 0.46,74 = 2) with f{Nm,X), we infer a baryonic mass density 
^Lya = 0.65Qb- Again, this should be considered a lower limit to Qi^ya subject to 
the uncertainties of Equation[T] More detailed comparisons of the opacity of the Lya 
forest with cosmological simulations reach a similar conclusion: the EUVB ioniza- 
tion rate exceeds that from quasars alone and Qi^ya — 0.95X2^ (e.g. l6ll l86l 18811 ). 
From a purely empirical standpoint, however, we must allow that as many as 20% 
of the baryons at z = 3 are still unaccounted for It is very unlikely that the difference 
lies in photoionized gas with A^hi < 1 c m~^ b ecause the mass contribution of the 
Lya forest peaks near A^Hi = lO'^^cm-2 UM- Instead, one would have to invoke 
yet another phase of gas, e.g. hot diffuse baryons. 

Ongoing surveys will fill in the current gap in f{Niii,X) for A^hi = lO^'^cm^^ 
down to at least A^hi = lO'^'^cm^^. This gas, which likely is the interface be- 
tween galaxies and the IGM, is very unlikely to contribute significantly to the bary- 
onic mass density. It may be critical, however, to the metal mass density at z ~ 3 



1 121 II 1511 . Furthermore, the partial Lyman limits (A^hi ~ 10 cm^ ) dominate the 
opacity of the universe to HI ionizing radiation and therefore set the "attenuation 
length" I42II . This quantity is necessary to convert the observed luminosity function 
of ionizing radiation sources into an estimate of Feuvb l54ll85ll . These observations 
may also provide insight into the processes of reionization at z > 6 and the escape 
fraction of star-forming galaxies. 



4.5 Hot Diffuse Baryons 

In the previous section, we demonstrated that the photoionized IGM may contain 
all the baryons in the z ~ 3 universe that are not associated with dense gas and 
stars. We also noted, however, that this conclusion hinges on the relatively uncertain 
intensity of the EUVB radiation field. As a consequence even as much as 30%Qh 
may be unaccounted for in the IGM. One possible reservoir in this category is a 
hot (r > lO^K) diffuse phase. Indeed, there are several processes already active at 
Z ~ 3 that can produce a reservoir of hot and diffuse gas. The hard, intense radiation 
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field from quasars and other energetic AGN processes (e.g. radio jets) may heat 
and ionize their surroundings out to several hundred kpc (e.g. [|20[l ). Similarly, the 
supernovae that follow star formation may also drive a hot diffuse medium into the 
galactic halo and perhaps out into the surrounding IGM (e.g. [2, 68]). Outflows of 
ionized gas are observed in the spectra of quasars and star-forming galaxies at these 
redshifts. On the other hand, gas accreting onto galactic halos is expected to shock 
heat to the virial temperature {T > lO^K), especially in halos with M > 1O"M0 
1 3411 . Altogether these processes may heat a large mass of gas in and around high z 
galaxies. 

Presently, there are very weak empirical constraints on the mass density and 
distribution of a hot, diffuse medium at z 3. Detecting this gas in emission is 
not technically feasible and there are very few observable absorption-line diagnos- 
tics, difficult to trace, even with absorption-line techniques. The most promising 
approac h use s the OVI doublet to prob e highly ionized gas associated with the IGM 
(e.g. ifsi I143II ) and galaxies |45, 144 1. These studies suggest a hotter component 
(T >~ lO^K) contributes on the order of a few percent Qi, at z~ 2. 

The mass density of an even hotter phase (~ 10^^^ K), meanwhile, is uncon- 
strained by observation. Ideally we would assess this component with soft X-ray 
absorption and/or emission but this will require a major new facility such as Con-X. 
We will take up these issues again in Section § l6.6l on the z ~ universe, where the 
warm/hot component is predicted to be substantial. 



Table 2 Empirical Summary of Baryons in the Universe 



Phase 


Temperature 


n{z = 3) 


Q 


(z = 0) 




(K) 


Location 


Estimate" 


Location 


Estimate" 


Stars 




Galaxies 


0.005 ±0.002 


Galaxies 


0.06 ±0.03 


Molecular Gas 


102 


Galaxies 


> 0.001 


Galaxies 


0.0029 ±0.0015 


Neutral Gas 


10^ 


Galaxies 


0.016 ±0.002 


Galaxies 


0.011 ±0.001 


Ionized Gas 


lO'* 


IGM 


>0.80 


IGM 


17+°- 

^•"-0.05 


Warm/Hot Gas 


10'' 


Galaxies? 


>0.01 


Filaments? 


?? 


Hot Gas 


10' 


?? 


?? 


Clusters 


0.027 ±0.009 



" Relative to an assumed Qi, value of 0.043. 



5 Baryons in the z ~ 1 Universe 

In contrast to the success that astronomers have had in probing baryons at z ~ 3, the 
knowledge of the distribution of baryons at z ~ 1 is quite limited. This is the result 
of several factors. One aspect is psychological: astronomers have pushed harder to 
study the universe at earlier times. This emphasis is shifting, however, with advances 
in IR spectroscopy (where the rest-frame optical diagnostics fall for z 1), the im- 
pact of the Spitzer mission (e.g. Dickinson's contribution), and the generation of 



Baryons: What, When and Where? 



19 



large spectroscopic samples of z ~ 1 galaxies (GDDS, DEEP, VIRMOS). The other 
aspect is technical. As with the z 3 universe, it is difficult or impossible with cur- 
rent facilities to build large samples of clusters to assess the host phase, or to survey 
the molecular phase with CO, or to measure HI gas with 21 -cm techniques. Fur- 
thermore, the majority of key absorption-line diagnostics (e.g. Lya, CIV; Table [T]i 
fall below the atmospheric cutoff and UV spectroscopy on space-based telescopes 
is necessary and expensive. Key emission-line diagnostics (Ha, [OIII]), meanwhile, 
shift into the near-lR making it more difficult to pursue large galaxy surveys. 

Surprisingly, the prospects are also limited for major advances in exploring 
baryons at z 1 during the next decade. The notable positive examples are (i) 
Sunyaev-Zeldovich experiments, which should establish the mass function of z ~ 1 
clusters, (ii) multi-object, near-lR spectrometers which will improve detections of 
star-forming galaxies and the luminosity function of early-type galaxies and (iii) the 
construction of mm telescopes (LMT, ALMA) that will enable searches for molec- 
ular gas. These advances will allow astronomers to characterize the dense, baryonic 
component at z ^ 1 . Similar to z 3 and z ^0, however, we expect that these com- 
ponents will be minor constituents. It is a shame, therefore, that there is currently 
no funded path to empirically constrain the properties of diffuse gas, neutral or ion- 
ized. Upcoming 'pathfinders' for the proposed Square Kilometer Array (SKA) will 
extend 2 1 -cm observations beyond the local universe but are unlikely to reach z = 1 . 
This science awaits an experiment on the scale of SKA. The classical absorption- 
line experiments must access hundreds of faint QSOsat z — 1 — 2 and so require a 
UV space mission probably on the scale of JWST 11371] . At present, however, there 
is no development funding path toward such a mission. In this respect, our igno- 
rance on the distribution and characteristics of the majority of baryons in the z 1 
universe may extend beyond the next decade. 



6 Baryons in the z ~ Universe 

In this section, we follow closely the footsteps of 14711. who presented a compre- 
hensive census of baryons in the local universe (see [46l| for an update). With this 
exercise in accounting, the authors stressed a startling problem; roughly 50% of the 
baryons in the local universe are missing! They further suggested that these missing 
baryons are likely hidden in a warm/hot (T « 10^ to lO^K) diffuse medium that 
precluded easy detection. These claims were soon supported by cosmological sim- 
ulations that predict a warm/hot intergalactic medium (WHIM) comprising 30 to 
60% of today's baryons IH, 31 ]. Motivated by these semi-empirical and theoretical 



studies, the search for missing baryons gained great attention and continues today. 
We review the current observational constraints on the distribution and phases of 
baryons in the local universe and stress paths toward future progress. As in previous 
sections, we roughly order the discussion from dense phases to the diffuse. 
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6.1 Stars 



Large-area imaging and spectroscopic surveys (e.g., 2MASS, SDSS, 2dF) have pre- 
cisely measured the luminosity function of galaxies in the local universe to faint 
levels H, 23]. The light emitted by galaxies in optical and IR bands, therefore, is 
well characterized. Establishing the mass budget of baryons in th ese st ars, however, 
requires translating the luminosity function into a mass function lll30ll . The typical 
procedure is to estimate a stellar mass-to-light ratio (M/L) for each galaxy based 
on stellar population modeling of the colors of the integrated light. These estimates 
are sensitive to the metallicity, star-formation history, and (most importantly) the 
initial mass function (IMF) of star formation. Because the stellar mass of a galaxy 
is dominated by lower mass stars even at early times, the analysis is most accurately 
performed in red passbands. 

The SDSS and 2dF galaxy surveys have led the way in this line of research; 
SDSS possesses the advantage of z-band imaging (A « 9000A) which reduces the 
uncertainty in the M/L estimate. 1670 analyzed the mass-to-light ratio for a large 
sample of SDSS galaxies, for which they report a luminosity function-weighted ratio 
< M/Lj >= 1 .5 with an uncertainty of w 20% and assuming a Kroupa IMF. Taking 
the |67] estimation of M/L with the SDSS luminosity function jH, |46] estimates 
a stellar mass density £2*{z = 0) = 0.0025, which includes an estimate of stellar 
remnants. Assuming a 50% uncertainty, we find £2*{z = 0) = (0.06 ±0.03)^2^,. 

Further progress in estimating the stellar contribution to Qi, is not limited by nei- 
ther statistical error nor depth of the luminosity function, but rather by systematic 
uncertainty, especially in the IMF. did not vary this assumption when estimating 
the error budget. Because sub-solar mass stars contribute most of the mass in stellar 
systems, uncertainties in the IMF can easily dominate the error budget. A conser- 
vative estimate of this uncertainty is « 50%, especially if we allow for a Salpeter 
IMF 

The subject of stellar masses and their sensitivity to IMF is worthy of a review 
article in its own right, and so it is well beyond the scope of our mandate. Here we 
only pause to note three recent developments. First, the uncertainty in total stellar 
masses is probably only a factor of 1 .5 - 2, even though the error may be larger for 
individual cases. Second, this topic is receiving more attention in stellar population 
synthesis models (e.g. Maraston, Bruzual and Chariot), which can now include IMF 
uncertainties in their budgets of systematic error. Finally, there have been recent, 
speculative suggestions that the IMF at high redshift may depart significantly from 
the Galactic form. At early times, generally hotter environm ents within galaxies 
and an elevated CMB may raise the typical stellar mass HEUEii. Though 



much theoretical work and experimental testing remains to be done in this area, it 
is possible that these ideas will ultimately lead to a more self-consistent picture of 
the evolution of stellar mass with redshift, with a consequent improvement in our 
understanding of the baryon budget. This is one of the key science goals of the 
JWST and ALMA facilities. These uncertainties aside, it is evident that stars are a 
sub-dominant baryonic component in our current universe. As an editorial aside, we 
note the irony that even the baryonic component that is easiest to see - the luminous 
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Stars - is difficult to weigh. This uncomfortable fact should be kept in mind later as 
we attempt to assess the unseen, hot diffuse phases of the IGM. 



6.2 Molecular Gas 



With the exception of the Galaxy 1113111 . the Magellanic Clouds 1116311 . and some 
nearby starburst galaxies [fsil]. there are very few direct detections of H2 in the local 
universe. Commonly these detections are of absorption in the FUV Lyman- Werner 
bands, and so are accessible only in optically thin, low-dust {Ay < 1) environments. 
H2 can also be detected in UV fluorescence lines from interstellar gas, but these are 
intrinsically very faint. Neither of these H2 reservoirs are likely to contribute much 
to the total cosmic mass budget of molecular gas (judging from the Milky Way). 
Because it is a homonuclear molecule with no dipole moment, H2 has only very 
weak quadrupole emission even if the density is high. The usual technique, then, is 
to track molecular gas with trace optically-thin species and infer the total molecular 
content through scaling factors calibrated against our Galaxy and others in the local 
universe. The most common indirect tracer is CO and its conversion factor is termed 
the "X-factor". In contrast with HI gas (see below), 'blind' surveys for CO have not 
yet been performed. Instead, researchers have estimated the molecular mass density 
by correlating CO mass with another galaxy observable (e.g. HI mass, infrared flux 
density) that has a well defined, volume-limited distribution function. T he convolu - 



tion of the two functions gives an estimate of the CO mass function 1I69LI175LI178I1 . 
Finally, this is converted to the H2 mass function with the X-factor which is gen- 
erally assumed to be independent of any galaxy properties. This latter assumption 
is known to be invalid for faint metal-poor galaxies and luminous starbursts, but 
variations are believed to be small for the L « galaxies that likely dominate the 
molecular mass density | [35| 



Taking X=A?(H2)//(C0)= 3 x 1020cm-2(Kkms)-i ifl?! ||69|l report X2h. = 



Phj = 2.2 X lO^/ivoMg, Mpc ^ for an IR sample of galaxies and Hi 7811 report Qa-, = 



1.2 X IO^/ivqMqMpc^^ for the optically selected BIMA sample. The factor of two 
difference in these estimates reflects the systematic error associated with sample 
selection. We adopt the mean value of the two estimates and assume a 30% uncer- 
tainty: i2H, = 0.0029 ±0.0009 f2t. The next major advance for surveying molecular 
gas in the local universe awaits a wide-field, deep, and blind survey of CO gas with 
depth and area comparable to current 21 -cm surveys. This is a difficult technical 
challenge and we know of no such survey that is currently planned. 



6.3 Neutral Hydrogen Gas 



Over the past decade, radio astronomers have performed bl ind surveys for 2 1 -cm 
emission to derive the HI mass function of galaxies (e.g. 11271 Il76l Il77ll ). The 



22 



Jason X. Prochaska and Jason Tumlinson 



modern version is the HIPA SS su rvey l87ll . an all Southern-sky survey for HI with 
the Parkes radio telescope. 117711 present the HI mass function from HIPASS and 
derive an integrated mass density X2hi = 3.8 x lO^^/iyg^ Including He, this implies 
a neutral gas mass density Qneut ~ (0.011 ± 0.001 )i2i. Although some uncertainty 
remains regarding the shape of the low-column ("faint-end") slope for the HI mass 
function, the total mass density is dominated by L w L* galaxies. Another important 
result of these surveys is the demonstration that very little HI mass, if any, associated 
with star-free systems. Ultimately, the ongoing ALFALFA survey 15211 will surpass 
the HIPASS sample in sensitivity and provide the definitive measure of the HI mass 
function at z ^ 0. In the future, facilities like the EVLA or Square Kilometer Array 
will be able to push such observations to z 1 or beyond. 

Interferometric observations can map the HI surface density profiles of low- 
redshift ga laxies and thereby infer the HI frequency distribution at A^hi > 10^'^ cm^^ 
I129I I179I1 . This distribution function is illustrated in Figure |7] Remarkably, it 
has nearly identical shape to the f{Niii,X) function derived for z ^ 3 galaxies 
from damped Lya observations (Figure|5]l. Both of these f{Niii,X) functions fol- 
low a roughly N^^^ power-law at lower column densities and exhibit a break at 
logA^Hi ~ 21.7. This suggests that HI gas is distributed in a similar manner within 
galaxies at z and z ~ 3. Perhaps even more remarkable is the fact that the in- 
tegrated mass density is also very comparable. The Q„eut value at z is at most 
three times smaller than the largest Q„eiit values derived from the z ~ 3 damped Lya 
syste ms (0.0045 vs. 0.0014) and actually coincides with the £2„ei,t value at z = 2.2 
1 1 14II . One might infer that X2„g„, has been constant for the past « lOGyr, which 
suggests that all of the_^as accreted into the neutral phase within galaxies is con- 
verted into stars 



or the gas 

(e.g. mi 



6.4 Intracluster Medium (ICM) 

The majority of galaxies in clusters are early-type with little or no neutral hydrogen 
gas. It may have been thought at first, therefore, that the baryonic content of clusters 
was almost entirely stars. X-ray observations, however, have revealed that clusters 
are filled with a hot plasma that extends to at least their virial radii (^ 1 Mpc). 
This intracluster medium (ICM) contains the baryons that have collapsed within the 
cluster structure and have been shock-heated to the virial temperature {T > lO^K) 
of the cluster Analysis of the X-ray surface brightness indicate that the ICM has 
an electron density profile «g(r) = no[l + (^'/^r)^]^^^^ with j3 w 0.6 and «o ranging 
from 10^' to 10^^ at rc of a few hundred kpc (e.g. 112211 1. In fact the ICM mass 
significantly exceeds the stellar mass inferred for the galaxies within the cluster 
From a cosmological standpoint, therefore, the ICM may represent a major baryonic 
component of our current universe. 

Detailed analysis of the X-ray observations show that the gas fraction, /g„s, or the 
ratio of the gas mass to the dynamical mass, is nearly constant with cluster mass, 
fgas = 0.1 1 ± 0.01 13J- Therefore, one can estimate the total baryonic mass density 
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of the ICM in all clusters simply by taking the product of fg^s with the cosmologi- 
cal mass density of clusters. A recent estimation of the latter comes from an X-ray 
selected sample of clusters with overlapping SDSS observations 112511 . Fitting these 
data to the 16211 dark-matter halo mass function for a WMAP5 cosmology 13711 . 
and integrating from high mass down to M„„„ = lO'^M©, we derive a mass density 
of clusters, ^cluster = 0.010. Adopting fg^s = 0.11, the ICM has a baryonic mass 
density Qjcm = 0.027i2f,. Uncertainties in this estimate are dominated by the un- 
certainty in Qciuster^ which are in turn driven by uncertainty in the power spectrum 
normalization specified by Cg. These errors translate into a roughly 30% uncertainty 
in Qciuster- This diffuse, hot component represents a non-negligible mass density in 
our modern universe. It raises the possibility that a similar medium associated with 
collapsed structures having M < M„n„ (i.e. groups and isolated galaxies) could con- 
tribute significantly to the mass density of baryons at z ^ I47i. i9Q1 . 
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f{Niii,X) distribution function of ionized Nm < 10^" cm ^, neutral Nui 



Fig. 

jq20-22^jjj-2^ and molecular Nm > 10^^ cm^^ gas at z ~ 0. The data at A^hi < lO' ^cm ^^ which 
correspond to HST/STIS observations of the low redshift Lya forest are taken from Il03 | ] . the dis- 
tribution function neutral, atomic hydrogen gas is taken from 21cm WSRT observations |l79[l , and 
the distribution function of molecular gas is estimated from BIMA CO observations Il78ll . 



24 



Jason X. Prochaska and Jason Tumlinson 



6,5 Photoionized Gas 



Summing the mass densities for the gas and stars of the preceding sub-sections, 
we estimate a total mass density of only « 0.10X2^. Similar to the high-z universe, 
therefore, we conclude that the majority of baryons aiz^O lie in diffuse gas outside 
of virialized structures. It is reasonable to speculate first that the photoionized com- 
ponent, i.e. the intergalactic medium, is also the major baryonic reservoir at z 0. 
As with the z ~ 3 universe, this component is best traced by Lya absorption using 
UV absorption-line spectroscopy. Even a visual inspection of low-z quasar spectra 
shows that the incidence of Lya absorption (.{X) is quite low relative to the z ~ 3 
universe. This sharp decline in is a natural consequence, however, of the ex- 
panding universe. Qua ntita tively, the decline in 1{X) is actually less than predicted 
from expansion alone il69ll because of a coincident decrease in the EUVB intensity 
which implies a higher neutral fraction and correspondingly higher incidence of HI 
absorption ifsill . 

The HI column density distribution function of the z ~ IGM has been char- 
acterized b y severa l groups using UV spectrometers onboard the Hubble Space 
Telescope 1331 110311 . In Figure |7] we show the Penton et al. f{N\ii,X) results for 
a sample with mean redshift z = 0.03. These results assume a Doppler parame- 
ter i>iGM = 25 km s^' for all absorption lines identified in their survey ( 12811 have 
obtained /jigm ~ 27 km s^' for their larger sample of 650 Lya absorbers). The re- 
sults are insensitive to this assumption for low column densities (A^hi < lO'^cm^^) 
because the Lya profiles are unsaturated but measurements at higher A^hi are sen- 
sitive to ^iGM and have much greater uncertainty. Similar to the z ~ 3 universe, the 
f{Niii,X) function is sufficiently steep that one expects the mass density is domi- 
nated by clouds with A^hi < lO'^cm^^. iEtIi have fitted the distribution function at 
low A^Hi with a power-law function /(A^hi,^) =^-^hi and report A = lo'''-^*' '' and 
/3 = — 1.73 ±0.04 for their recent sample with (z) = 0.14, and there is significant 
degeneracy between the two parameters. 

We may estimate the IGM bary onic mass density at z using the formalism 
described in § 14.41 following lll34ll . In addition to f{Niii,X), this analysis requires 
an estimate of the HI photoionization rate F and the gas temperature Tigm- The de- 
pendences of the mass density on these quantities are weak {Q^ya °= F'/^T^/^) but 
so too are t he em pirical constraints. Estimates of F range from 0.3 to 3 x 10^'^ s^' 
lfl36[ [l39[ [iToll and we adopt F = 10 s ' based on calculations of the com- 



posite quasar and galaxy background estimated by Haardt & Madau (CUBA; in 
prep). The characteristic temperature of the IGM is constrained empirically only 
through the observed distribution of line widths, ^igm- Several studies report a char- 
acteristic line width at z of biGM ~ 25 km s^' which implies Tigm < 38,000K 
I33l[l0fll 10211 . This is an upper limit because it presumes the line widths are ther- 
mally and not turbulently, broadened. Indeed, 13311 have argued from analysis of 
numerical simulations that the observed /jigm distribution is in fact dominated by 
non-thermal motions. These same simulations suggest a density-temperature rela- 
tion TiGM ~ 5000(p /p)''^K (the so-called "IGM equation of state"). 
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Combining F = 10 s ' and the | 33n p — T "equation-of-state" with the il03 1 
estimation for f{Nm,X) and the 113411 formalism, we estimate the baryonic mass 
density of the z - Lya forest (A^hi = lO'^-^^^^'^cm^^) be X2Lya = 0.0075 = 
0.1712/,. From their sample of 650 Lya absorbers at (z) = 0.14, |28J obtain Qiya = 
0.0131 ±0.0017 = 0.29X2^,. This range is still several times smaller than the mass 
density estimate for the IGM at z ~ 3 (Figure|6]l. The immediate implication is that 
the IGM may not be the majority reservoir of baryons at z ~ 0. Indeed, the missing 
baryons problem hinges on this result. It is crucial, therefore, to critically assess the 
uncertainty in this estimate, in particular whether current observational constraints 
allow for a significantly larger value. 

First, consider the photoionization rate F which is constrained to no better than 
a factor of a few (at any redshift). As noted above, Qi^ya is largely insensitive to 
the photoionization rate (Qhya °^ F^/^). Furthermore, our adopted value lies toward 
the upper end of current estimations which leads us to conclude that a more precise 
measure of F will not markedly increase the baryonic mass density inferred for the 
IGM at z - 0. 

Another consideration is the temperature of the IGM, although again the mass 
density is not especially sensitive to this quantity. Even if the Lya forest were uni- 
formly at TiGM = 2 X 10'*K, we would derive only a 50% higher baryonic mass 
density. It is evident that modifications to F and Tigm will not make a qualitative 
difference in Q^ya^ though these effects in combination could increase it by a fac- 
tor of 2. Turning to the HI frequency distribution, Iil03il and |[33l emphasize that 
f{Niii,X) shows no break from a power-law down to the sensitivity Umit of the 
HST/STIS observations, Nii,„ — lO'^'^cm^^. It is reasonable, therefore, to extend 
the distribution function to lower column densities. Taking A^/,„, = 10"^cm^^ with 
the power-law exponent fixed at j3 = — 1 .65, leads to a 2 x increase in £2i,ya- Finally, 
both cosmological simulations and some observational analysis suggest a steeper 
power- law (j3 « -2) ||33, 82, 100] which, in turn, may lead to an enhancement in 
^Lya- If we allow for a combination of these effects, it is possible to recover Qi^ya 
values exceeding 0.5X2^,. For example, the following parameter set - F = 10^'^ s^', 
TiGM = 2 X lO'^K, J3 = -1.8, Nii,n = 10^'-^ cm^^ _ yields Qi^ya = 0.6X2^,. 

On empirical grounds alone, then, the current observations do not rule out a sce- 
nario where the photoionized Lya forest contains most of the baryons not locked 
in denser phases, which we estimate at ~ 0.1X2^. This would require, however, a 
warmer IGM than predicted by numerical simulations and a f{Niii,X) function that 
remains steep to A^hi ~ lO" cm^^. The latter point may be tested with new observa- 
tions of the IGM using the planned HST/COS instrument, especially if it achieves 
S/N > 50 per pixel. Given that the uncertainties must be pressed in concert to their 
limits for the budget to be complete, it is also reasonable to conclude that the IGM 
does not contain all of the missing baryons at z ~ 0. In the next section we will take 
up the question of what could fill this gap, assuming it exists. 
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Fig. 8 Estimates of the baryonic mass density (relative to £2^) for various phases of baryons in 
the c ~ universe. On the left side of the figure we show the mass densities estimated from the 
photoionized Lya forest and the neutral gas in low z galaxies. One notes that the sum of the central 
values for the various components is significantly less than unity. This suggests that a significant 
mass component is missing from this census, e.g. the warm/hot intergalactic medium (WHIM). 
Even adopting the maximum values for each component, one may require an additional baryonic 
phase to contribute at z = 0. 



6.6 The WHIM 

In Figure |8] we present the z ^ baryonic mass density estimates for the phases 
described in the previous sub-section. A simple summation of the central value es- 
timate of each componenj^ yields a total Qstms + ^Hi + ^neut + ^iCM + ^Lya = 
0.21 Qj,- Therefore, having followed the footsteps of FHP98, we reach a similar con- 
clusion: current observations have not revealed roughly 70% of the baryons in the 
Z ~ universe. The observational techniques described in the previous section are 
sensitive to baryons in luminous matter, hot {T > lO^K) and diffuse gas, molecular 
and atomic neutral gas, and diffuse photoionized gas. From a practical standpoint, 
this leaves one obvious phase: warm/hot (10^ K<T < 10^ K), diffuse gas. 

Indeed FHP98 identified this phase as the likely reservoir for the remainder of 
baryons and suggested it would be associated with galaxy groups and galactic halos 
(see also ||90ll). Cosmological simulations of the low redshift universe, meanwhile, 



^ Although there lies a gap in our knowledge or gas with A'hi ~ 10'^ cm ^, it is unlikely that this 
gas contributes significantly to Qf,. 
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argue that a warm/hot phase exists yet external to galaxies, i.e. within the intergalac- 
tic medium lllSlblll . This warm/hot intergalactic medium (WHIM), believed to be 
heated by shocks during large-scale gravitational collapse, is predicted by the sim- 
ulations to comprise roughly 50% of baryons in the present universe. Bregman ll3 l 
has recently reviewed the observational and theoretical evidence for the WHIM at 
Z^Q. We refer the reader to that manuscript for a broader discussion. We highlight 
here only a few points with emphasis on absorption-line observations. 

The WHIM gas, as conceived in numerical simulations or as gas reservoirs of 
galaxy groups, has too low density and temperature to be detected in emission by 
current X-ray instrumentation and even the next generation of X-ray facilities would 
be unlikely to detect the most diffuse gas. The neutral fraction, meanwhile, is far 
too low to permit detections via 21cm observations. For these reasons, searches for 
WHIM gas have essentially been limited to absorption-line techniques. The current 
approaches include (1) a search for the wisps of HI gas associated with the WHIM 
and (2) surveys for ions of O (and Ne) which trace warm/hot gas. 

The first (and most popular) approach has been to survey intergalactic along 
quasar sightlines. This is primarily driven by observational efficiency; oxygen ex- 
hibits a strong OVI doublet at ultraviolet wavelengths (1031, 1037 A) that can be 
surveyed with high resolution, moderate signal-to-noise observations. To date, these 
have been acquired by spectrom eters on the HST and FUSE satellites. In collisional 
ionization equilibrium (CIE; 1 53 , 150ll ). the ion is abundant in gas with T = 10^ 
to 10*'K and therefore may trace the cooler WHIM. A hard radiation field, how- 
ever, may also produce substantial ions in a diffuse gas. Therefore, the de- 
tection of OVI is not definitive proof of WHIM gas. Tripp and collaborators were 
the fir st to demonstrate that intergalactic OVI may be detected along quasar sight- 
lines 1 156l 157 1. Surveys along tens of si ghtlines h ave now been comprised with 
HST/STIS and FUSE observations ll^lHflillisl which reveal a line density i{z) 
of absorbers £{z) w 10 per unit redshift path down to Wio3i > 50mA and £{z) « 40 
down to Wi 031 > 10mA. 

Intriguingly this incidence and the observed equivalent width distribution are in 
broad agreement with the predictions from z ~ cosmological simulations lll7ll . At 
the least, this should be considered circumstantial evidence for WHIM gas. On the 
other hand, at least 30% of OVI systems s howing OVI absorption also exhibit co- 
incident CIII absorption 128 . miiilEsl. The 0+5 and C++ ions cannot coexist 
in CIE; their coincident d etecti on either indicates a multi-phase medium or a pho- 
toionized gas. Tripp et al. 115811 have also emphasized that « 40% of OVI absorbers 
show HI absorption with nearly identical line profile (Figure|9l see also 115311 ). In 
their analysis, Danforth & Shull fl^ have placed this fraction at < 10%, perhaps 
owing to their choice to not divide the absorbers into velocity subcomponents. Thus 
there is agreement that at least some, if not all, of these absorbers have temperatures 
ofr<10^K, i.e. photoionized materia^. In other cases, the OVI absorption is offset 
in velocity from the strongest HI absorption and there is little diagnostic constraint 
on the gas temperature. In only a few cases (e.g. the z = 0.1212 absorber toward 



* Although one could invoke undetected, broad Lya absorption in these cases, this would not 
predict the observed identical line-profiles for OVI and Lya absorption. 
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Fig. 9 Plot of the Lya, Lyj8 , and OVI line-profiles for an OVI absorber at z = 0. 12 toward 3C273 
(taken from Il58ll ). Note the high similarity between the OVI and HI profiles. This indicates the 
gas is not thermally broadened and requires a gas temperature T < lO^K. In turn, this implies 
the material is predominantly photoionized, not coUisionally ionized and is therefore unlikely to 
correspond to warmyhot intergalactic gas. Roughly 40% of OVI absorbers share these properties. 
It remains an open question as to whether the majority of observed OVI gas corresponds to the 
WHIM. 



HS1821+643; 115511 ) does one detect broad Lya absorption aligned with the OVI 
absorption that clearly favors coUisionally ionized gas with T > lO^K. 

In support of the view that the OVI absorbers are predominantly hot and coUi- 
sionally ionized, Danforth & ShuU 112 811 offer the arguments that (1) the column- 
density and linewidth distributions of the two high ions OVI and NV differ signif- 
icantly from those of commonly photoionized species, such as CIII/IV and HI, and 
(2) the observed systems imply clouds with an unrealistically wide range of ion- 
ization parameters in close proximity to yield OVI and lower ions. Of the known 
OVI systems, the majority have not been confirmed by multiple investigators as 
truly arising from collisional ionization, so there is presently no consensus that OVI 
in fact measures the r > 10^ WHIM. The upshot of these ongoing controversies is 
that while OVI can be a tracer of WHIM, its origins are ambiguous and further work, 
including detailed cross-checking by the interested groups, is needed to assess the 
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fraction of OVI absorbers that are truly hot and colhsionally ionized, and from there 
to calculate their contributions to the baryon budget. 

It may in fact be premature to identify this OVI-traced gas with the intergalactic 
medium. Although the incidence of OVI absorbers is too high for all of them to be 
associated with the halos of L « L* galaxies, there are sufficient numbers of dwarf 
galaxies to locate the gas withi n their hal os or "zones of metal enrichment" if these 
extend to w 150-250kpc lEl [iH [160I1 . Indeed, searches for galaxies Imked with 
OVI absorbers have indicated a range of as sociations frorn indiv idual galactic halos 



to galaxy groups to intergalactic gas 11251 ll 17L Il48l |I54 II62I1 . We await a larger 



statistical sample and detailed comparison to numerical simulations to resolve this 
issue l48ll . The upcoming Cosmic Origins Spectrograph on HST will make a large 
and important contribution to solving this problem. 

Setting aside the controversial aspects of OVI, one may speculate on the mass 
density of the OVI-bearing gas i2ovi using similar techniques as for the Lya forest. 
Because OVI is simply a tracer of the material, one must assume both the metallicity 
and an ionization correction for the gas. Current estimations give £2oyi ~ 0. 1X2^,, as- 
suming that all the OVI detections arise in collisionally ionized gas |28]. Therefore, 
even in the case that the OVI-bearing gas traces only WHIM material, it is evident 
that it comprises a relatively small fraction of the baryon census. 

Ideally, one would assess the mass density and spatial distribution of WHIM 
gas through observations of hydrogen which dominates the mass. This is, how- 
ever, a difficult observational challenge. Assuming CIE, a WHIM 'absorber' with 
(an optimistic) total hydrogen column of Nh = lO^'^cm^^ and T = lO^K will 
have an HI column density of A^hi = lO'^^ ''^cm^^ and Doppler parameter bni = 
129 km s-i(r/10^K)i/2. This gives a peak optical depth Tq = 1.5 x IQ-^iNlLyafhya) 
0.25 and a total equivalent width Wcya — 137mA. The detection of such a line- 
profile requires at least moderate resolution UV spectra at relatively high S/N and 
also a precise knowledge of the intrinsic continuum of the background source. 

Detections of line profiles consistent with broad Lya absorption have now been 
reported in the few HST/STIS echelle spectra with high S/N lilH^ . Because 
these detections lie near the sensitivity limit of the spectra, it is not possible to 
exclude the possibility that these broad features arise from the blends of several 
weak, narrow absorbers. Furthermore, independent analysis of these same datasets 
have questioned most of the putative detections (Stocke, priv. comm.; Danforth & 
Shull, in preparation). These issues can only be resolved by higher S/N observa- 
tions and, ideally, analysis of the corresponding Ly)3 profiles. We look forward to 
the installation of HST/COS which should enable such observations. Even if one can 
conclusively demonstrate the presence of warm/hot gas with broad Lya absorbers, 
there are limitations to directly assessing the mass density of this gas from empiri- 
cal observation. First, studies of numerical simulations reveal that broad Lya lines 
arise from gas whose kinematics are a mixture of thermally and turbulently broad- 
ened motio ns m aking it difficult to derive the temperature and obtain an ionization 
correction 112311 . Second, these same models predict that photoionization processes 
are significant for some gas with T < lO^K and this contribution to the ionization 
correction cannot be assessed from HI observations alone. Finally, the detection of 
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broad Lya features becomes prohibitive for gas with T > 10''^K (i.e. Tq < 0.05 for 
Nh = 10^^ cm^^), which may account for the majority of WHIM material by mass 
(e.g. 13111 ). In short, one is driven to statistical comparisons between numerical sim- 
ulations and observations of broad Lya absorption and, ultimately, may be limited 
to a minor fraction of the WHIM. 

At r > lO^^K, the ion fraction of is negligible and broad Lya lines are 
generally too weak to detect with previous or planned UV spectrometers. To more 
efficiently trace this gas, we must turn to higher ionization states of oxygen and 
other elements. One option is the NeVIII doublet whic h pro bes gas to T w lO^K 
and is accessible with UV spectrometers for z > 0.25 lll32ll . The most sensitive 
tracers of T > lO^K gas, however, are resonance lines of OVII (21.602 A), OVIII 
(18.969 A), and the NeVIII doublet (770,780 A). While NeVIII may be detectable 
at z > 0.5 with COS, for current X-ray satellites the oxygen lines are a demanding 
observational challenge; one is limited by both spectral resolution and instrument 
collecting area. Nicastro and collaborators have addressed the latter issue (in part) 
by monitoring X-ray variable sources that show occasional, bright flares. During a 
flare event, one obta ins ta rget-of-opportunity observations to yield a relatively high 
S/N dataset 11411194 117211 . While this program has been successful observationally, 
the p urported detections of intergalactic OVII and OVIII remain controversial i65l 
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While we are optimistic that these programs will reveal a few indisputable 
detections before the demise of the Chandra and XMM observatories, it is evident 
that a statistical survey for OVII and OVIII awaits future instrumentation (e.g. Con- 
X, XEUS). Because OVI may indicate cooling gas within a hot medium (or hot- 
cold int erface s), the numerous detections of OVI can serve as signposts to the 'true' 
WHIM lfl40ll . 

Ultimately, however, the challenge will remain to convert these tracers of hot 
gas into constraints on the mass density and spatial distribution of WHIM gas. At a 
minimum, this requires assumptions about the metal-enrichment of the intergalactic 
medium. 

In summary, we find the current empirical evidence for the WHIM to be intrigu- 
ing but far from definitive. One must consider additional means for exploring the 
WHIM that would be enabled by new, proposed or planned instrumentation. One of 
these is to explore the line emission (e.g. Lya, OVI, OVII) from the denser WHIM 
gas ||7l[l3], which could be accomplished with UV and X-ray missions currently in 
design (e.g. IGM,XENIA). If combined with absorption-line studies, one may map 
out the mass density and spatial distribution of this gas at low z. Finally, it may 
be possible to cross-correlate the Sunyaev-Zeldovich (SZ) signal of the CMB with 
WHIM structures to infer the density and spatial distribution of this ionized gas 115 511 . 
This might be accomplished with currently planned SZ experiments. 
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Fig. 10 For specific missions or categories of facility, the colors represent the role they are ex- 
pected to play in measuring the budget of baryons in each phase discussed in the text. 

7 Future Facilities and the Baryon Census 

In the spirit of the meeting "Astrophysics in the Next Decade: JWST and Concur- 
rent Facilities", we consider how the baryon census will be advanced by the facili- 
ties in development now for operations in the next decade. To stimulate discussion 
at the conference, we composed a chart showing the applicability of each mission 
or facility to the important reservoirs of baryonic matter. The list of missions un- 
der consideration was obtained from the conference website. This chart appears in 
Figure [TOl with the contributions of each mission color coded for major (green), 
supporting (yellow), minor (red), or no (grey) role. These assignments were made 
in the idiosyncratic judgment of one of us (J. T.) and should not be taken too liter- 
ally. Each mission has its own, presumably compelling, scientific justification that 
in most cases does not explicitly include the baryon census (though JWST and Con- 
X are notable exceptions). The assignments were based primarily on wavelength 
coverage, spectral and spectroscopic resolution, field of view, and a review of the 
mission's public science case where available. 

The overlap between missions and baryonic components raise some interesting 
points. First, the future looks relatively bright for the baryon census as a whole. 
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JWST, Con-X, and ALMA are the three most important missions for resolving un- 
certainties in the census of baryons in stars, hot gas, and cold gas respectively. These 
are the areas of greatest uncertainty today, so this is somewhat encouraging for the 
future. However, only two of these facilities are in development phases: Con-X is 
still under study and could not reasonably be expected to fly until late in the next 
decade or after 2020. 

The "Assembly of Galaxies" is one of the four top-level science drivers of JWST, 
which will have the sensitivity and wavelength coverage to improve on stellar mass 
estimates at high redshift and the ability to measure IMFs by star counts over much 
larger local volumes than has been possible to date. As this is the most uncertain 
step in the derivation of the cosmic stellar mass budget, assessing the variation of 
the IMF from place to place and over time should be a primary goal of galaxy 
formation and evolution studies with JWST. 

It looks as though the next decade will bring major advances in the study of the 
cold neutral and molecular phases of baryons over a wide redshift range, thanks to 
ALMA and eVLA. As mentioned above, the greatest need in this area is for a more 
thorough testing of the conversion factor from CO mass to H2 mass in a broad range 
of astrophysical environments. 

One pleasant surprise from this exercise is the notable versatility of large, ground- 
based optical telescopes (Keck and VLT today, and the 20 — 30 m ELTs in the fu- 
ture). Though they can lead the way for direct studies of stellar masses over a wide 
range of redshifts, they also support investigations whose primary use is in another 
wavelength range. Often this role requires them to obtain supporting data such as 
spectroscopic redshifts for high-z galaxies, optical counterparts for X-ray or radio 
sources, and the measurement of dynamical masses by high-resolution studies of 
linewidths in absorption or emission. It is likely that this versatility reflects the long 
history of optical astronomy, where the most mature technologies and techniques 
are available, and their widespread availability compared with the scarce resource 
of space-based observing time. 

This exercise also highlights a broad trend we may expect from astrophysics in 
the next decade. Since the baryon census encompasses many diverse astrophysical 
environments, from stars to the diffuse IGM, and requires a wide range of observa- 
tional techniques, the capability of a facility to address the baryon census reflects its 
ability to address astrophysical problems in general. With this in mind, we note that 
many of the space missions in question are tailored to specific scientific questions 
rather than a broad, community-generated scientific program. While it is difficult 
to see how a mission like LISA could have collateral uses beyond the detection 
of gravity waves from a few astrophysical sources, even missions carrying more 
conventional instruments like JDEM, GAIA, and SIM are fine-tuned to specific sci- 
entific cases. This trend toward specialization, even among flagship-scale missions, 
will evidently be a hallmark of astrophysics in the next decade. 

Finally, we note that perhaps the greatest need for new capabilities to address the 
"missing baryons" lies in the X-ray, and that there is no realistic chance for Con-X 
or something like it to launch before the very end of the next decade. Thus it is 
possible that astronomers will, by 2015, see the "first galaxies" (JWST), measure 
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DE equation of state (LSST+JDEM), weigh molecules at z 6 (ALMA), discover 
gravitational waves (LISA), find earthlike extrasolar planets (SIM+GAIA), and still 
not know where all the ordinary matter is or what phase it is in. This long interval 
would place the solution to the missing baryons problem at least 20 years away from 
when it was first recognized. This may be one of those cases where the urgency of 
performing an observation is tempered by the conventional wisdom that we already 
know the answer (the WHIM). But theory has been wrong before, and the only 
way to know is to do the experiment. If the last decade of astrophysics is a reliable 
guide, the next ten years will bring a few surprises. The what, when, and where of 
the baryons may well be among them. 
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